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SYMBOLS 


f frequency, Hz 

j imaginary operator, 

k ratio of sample period to frequency analyzer step period 

s Laplace variable 

T computer step size and sample period, sec 

Tjq time lag of input/output operation, sec 
X(s) Laplace transform of continuous variable 
X*(s) Laplace transform of sampled variable 

T time constant of first-order system, sec“^ 

^ phase, radian 

difference between predicted and observed phase, deg 
w frequency = 2Ttf, rad/sec 

Subscripts 

FA frequency analyzer 

FI low pass filter 

F2 pure time lag "filter" 

DAC digi tal-to-analog converter 

Fre analog prefilter of ADC 
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A TIME LAG STUDY OF THE VERTICAL MOTION SIMULATOR COMPUTER SYSTEM 


William B. Cleveland 
Ames Research Center 


SUMMARY 


A study was performed to evaluate an experimental method to determine 
time lags in real-time computer systems as the one associated with the Vertical 
Notion Simulator at Ames Research Center. The approach was to use an ordinary 
frequency analyzer to measure the phase difference between inputs and outputs 
of the computer system. !1* various elements of the program and computational 
architecture were modele.. Various factors, such as computer frame time and 
input frequency, were varied so that they were representative of the opera- 
tional use of the simulator facilities. Experimentally determined results 
were compared with predictions derived from the simulation models. The results 
indicate that the frequency analyzer c^n be readli;, used to evaluate time lags 
in systems of this type. Differences between predicted and measured phase 
values Indicate that the hardware and software imparts a time lag of about 
.S msec to this facility. 


INTRODUCTION 


The effect of the computer system on the results of simulations has many 
facets and has been the subject of many investigations. The results can be 
affected bv the quality of the modeling of the system, by the method of dis- 
crete implementation of continuous models, by the quality of the progr.inuning 
technique, and by the hardware configuration of the computer system togetl’.et 
with its Input/output (10) system. This Investigation addresses itself to an 
analysis of the computer and 10 system associated with tlie Vertical Motion 
Simulator at Amt's Research Center. This simulator, because of its large ver- 
tical travel of *10 m, is attractive to researchers and program managers 
interested in the performance of helicopters and modern vertical and short 
takeoff and landing (V/STOL) aircraft. Simulation of these veliicles is more 
difficult than that of other aircraft, especially transport aircraft, hi'cause 
of the size of the aircraft model and, mc^re importantly, tlie higher natural 
frequencies of both the aircraft and their associated control systems. The 
combination of size and frequency requires that a large amount of computation 
be performed within a short time period. Therefore, it becomes important 
rid the simulation of unnecessary time lags since time lags ha ’e a destabiliz- 
ing effect on simulations of aircraft often already burdened with stability 
prob lems . 

Tile VMS computer system consists of two computers and an 10 system. 
Because of the number and location of the flight simulators (and computer 
facilities) at Ames, a policy of associating a small computer with the simu- 
lator is evolving. The "iiost" or large computer that solves the aircraft's 


(a) What Is the minimum turnaround time In the system? In other words. 
If the ADC samples a signal that is then used in the digital computer somehow 
and then is output via a DAC, what is the minimum time to perform this? 

(b) Are there any unknown time lags in the system? When one models the 
system and the contents of digital simulation, one has the means to predict 
the time lags. Do theoretical and actual results agree? 

(c) Are there simple means to measure the time lags in the digital simu- 
lation? The answer to this question, in general, is "no" if one is Interested 
in measuring the time lags of the blocks such as those in figure 1. These 
timings are merely incidental information for the systems analysis of Interest 
anvway. However, one can use a frequency analyzer to measure the total phase 
difference between the ADC's input anJ the DAC's output. This phase can be 
interpreted as total system time lag. 

These questions form the background for the general objective of the 
investigation which was to determine if a frequency analyzer can be used 
experimentally to provide verification or validation for a simulation computer 
conf igurat ion. 


Approach 

The method of satisfying the objectives consisted of comparing experi- 
mi*ntal and theoretical results. Figure 2 Illustrates the major elements for 
which linear models were derived. In the case of the "computer program," 
simple linear systems were selected for ease of operation as well as for 
clarity of interpreting results; in particular, a one-T delay was employi-d 
as one of the cases. 



GAIN PHASE 


Figure 2.- Frequency analyzer technique for determination 

of system time lags. 
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DESCRIPTION OF FACILITY 


The facility of interest is the computer and 10 aaaoclated with the 
Vertical Motion Simulator (VMS) at Ames. The capabilities of the simulator 
are described in reference 2. Within this section, the computer and 10 hard- 
ware and software are described to the extent that the models needed for the 
analysis may be derived. 


Computer and 10 Operation 

The relevant aspects of the computers and 10 systems are the method of 
timing and control of the various processes. The ADCs of the 10 section run 
asynchronously while the two computers operate synchronously. While the 
details of the 10 are given in reference 3, it is useful to outline its opera- 
tion here. The 10 system consists of Individual circuit boards for the ADC 
function (as well as for the DACs and discrete functions). Each board has a 
microcomputer to control 16 channeln of ADC. The controller and ADCs are all 
self-contained on the one board. If more than 16 ADCs are needed for the 
simulation, they arc added board by board to the 10 system In increments of 
16 ADCs. 

The microcomputer runs an ADC and a 16-channel multiplexer to effectively 
obtain 16 channels zi ADC. The operation consists of the microcomputer run- 
ning the function at 50 usee per ADC conversion or 0.8 msec for the 16 ADC 
channels. The results of the conversion are stored in a buffer on the circuit 
card. The buffer Is updated or has fresh information every 0.8 msec. Since 
each ADC board operates independently, each set of 16 ADCs Is refreshed at 
this rate also. Consequently, no matter how many blocks of ADCs are used, 
tliere is no more chan 0.8 msec of time skew between any two channels. 

Timing of Che overall computer system is controlled by a clock associated 
with the PDP computer. Figure 3 illustrates the process in which the ADCs, 
DACs, and Che host computer are coordinated by the PDP. Several important 
points are to be considered. The first is the order of the ADC-DAC-compute 
schedule. The implication is that, in a given period, Che DAC cannot output 
data that arc a function of the ADC input from that same period. /\n alternate 
schedule sometimes used in process control applications is the ADC-compuCc- 
DAC-compuCe schedule. The control laws that contain lead are calculated in 
the first computer period and are output via the DACs as soon as possible to 
minimize phase lag. Any remaining controls are calculated in the second com- 
puter period for an operation much the same as the ADC-DAC-compute schedu. e. 
This latter sequence of 10 and compute is not used in the simulation system 
at Ames. 

Another observation is that the individual timing associated with sending 
data from the ADC to the PDP, from the PDP to the host, etc., is not important 
from a systems point of view. So long as the input is sampled by the ADC at 
the first clock pulse and the output is converted by the DAC at the subsequent 
clock pulse, what went on in the computer system to make it all happen is not 
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1. ADC VALUES ARE OBTAINED WITHIN 0.8 msM IN ADVANCE 
OF THE CLOCK PULSE. THE LOCAL COMPUTER OBTAINS 
THE ADC DATA FROM THE i/O SYSTEM AND SENDS IT TO 
THE HOST COMPUTER AS INPUTS TO THE SIMULATION 
EOUATIONS. 

2. THE LOCAL COMPUTER LOADS AND INITIATES THE DAC 
CONVERSION IN THE I/O SYSTEM. 

3. THE HOST COMPUTER SOLVES THE EQUATIONS THEN 
OUTPUTS TO THE LOCAL COMPUTER WHICH READIES 
THE DAC DATA FOR THE NEXT CYCLE. 

Figure 3.- Timing and functional description of control cycle. 


important. Naturally, minimizing the "overhead" time required by all the sys 
terns to operate is important so that the simulation computing time will be 
sufficient; it is advantageous to minimize T to maximize accuracy of the 
simulation. 


System Models 

The system may be broken into the four components shown in figure 2, 
that is, the ADCs, the computer, the DACs, and the frequency analyzer. The 
model of each component contributes to the understanding of the system. 

ADC <vodel- The ADC model has two parts of Interest: the ADC prefllit r 

and the ADC Itself. The prefilter was included in the VMS 10 system to i . . 
out bO-Hz noise; its corner frequency is set to 20 Hz for this purpose. ;'ht- 
filter may be switched in or out of the system as desired. The phase of th> 
filter for input frequency ui below the corner frequencv is 

ifpre ■ -u)/40ti , u < wq 

From this, the equivalent pure time lag may be calculated as 
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important. Naturally, minimizing the "overhead" time required by all the ays 
terns to operate is Important so that the simulation computing time will be 
sufficient; it is advantageous to minimize T to maximize accuracy of the 
simulation. 


System Models 

The system may be broken into the four components shown in figure 2, 
tliat is, the ADCs, the computer, the DACs, and the frequency analyzer. Ttie 
model of each component contributes to the understanding of the system. 

f\DC model- The ADC model has two parts of Interest: the ADC preflltir 

and the ADC Itself. The prefilter was Included in the VMS 10 system to f'. . 
out bO-Hz noise; its corner frequency is set to 20 Hz for this purpose. ITu- 
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From this, the equivalent pure time lag may be calculated as 
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Tpre ■ “ l/40ir ■ 0.008 sec 

Thus filtering the high frequencies on the input results in a pure time lag 
of 8 msec. 

The ADC converter model shown in figure 4 is made up of time delay plus 
a sampler operating at the rate of one sample for every .-eriod of T seconds. 
The sampler is the standard representation of an ADC, but the addition of the 
time lag of Tjq seconds is not. The lag is Introduced for two reasons. One 
reason is that the 0.8-msec period for sampling 16 ADCs is extended by the 
actual data transfer process from the ADC microcomputer to the PDP computer. 

If the sampling rate T is 40 msec, 50 conversions of the 16 ADCs have 
occurred every time the PDP computer "samples" the ADC unit. Since the PDP 
may sample at any point in the ADC's operation, it is possible that the con- 
verted ADC data in the microcomputer may be as much as 0.8 msec old, exactly 
in synchronization, or somewhere in between. An average time lag of 0.4 msec 
In the ADC is a reasonable assumption. 


PREFILTER 



Figure 4.- ADC system model containing a tourth-order 
Bessel filter and conversion model. 


The second portion of the 10 time lag is due to the time it takes the PDP 
to receive the ADC data from the RIOU and send it on to the Sl>’ma computer in 
addition to the output data transfer from the PDP to the RIOU. This time has 
been estimated to be on the order of 3 to 4 msec but no timing has been per- 
formed. Since timing data are not available, a pure time delay of Tjq is 
assumed to be present with its magnitude to be determined from the differences 
between experimental and theoretical timing. 

The model of the sampler itself may be obtained from reference 4: 

00 

X*(s) ‘ I 21 - Jn2Ti/T) 

n»® 

In this experiment, the repeating spectra generated by the sampler have no 
effect on the experimental results. Although they are the source of the 
"staiicase" effect in a DAC's output, it is shown in the appendix that, given 
a frequency analyzer that operates on the basis of forming Fourier series 
coefficients, the repeating spectra have no effect on the coefficients when 
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Che ratio of the sampling frequency Co Che input frequency being analyzed 
divides evenly into 2C0. This condition was observed in Che experimental 
setup. Thus Che ADC model may be simplified to 

X*(s) - ^ X(s) 

Computer model- From the considerations discussed in the section on com- 
puter and 10 operation, no model is required of Che computer hardware Itself. 
What is needed is the model of the computations performed. Two different 
dynamical models were simulated for the test ~ one a first-order system of 
unity gain and time constant T, and the second, a pure time lag of T seconds. 

The first-order system was mechanized using the zero-order hold method 
resulting in the theoretical phase: 

_i -sin iijT 

FI cos oj, - exp(T/T) 

The pure lag system produces the phase: 



Both filters were run under the Ames software (described in i'efs. 4 and 
5), In which portions of the calculations may run at speeds that are fractions 
of the basic or nominal speed. For this test, both filters were simulated 
using periods T and 2T. The effect of operating the software at the slower 
speed, equivalent to a period of 2T, is to place a sampler in the "slow" path 
which, In effect, performs the function of an ADC running at period 2T. The 
computations produce a new result every 2T seconds which effectively makes 
the DACs operate at a period of 2T also. 

Wlien the experimental tests are conducted, representations of both fil- 
ters art operating at both fast and slow speeds; thus, there are four computer 
model cases. 

DA' model- The dlgitai-to-analog converters are conventional zero-order 
hold devices. Like the ADCs, the DACs are packed 16 to a circuit board. Fig- 
ure 3 illustrates the timing in which the PDF computer loads buffers of the 
DAC circuit and initiates the actual conversion. Since the 10 lag was dis- 
cussed in the ADC section, the DAC operation may be modeled in the conventional 
manner. The phase relationship for the zero-order hold DAC (taken from iii. 6) 
is 


*DAC ■ -wT/2 

This phase is Interesting since It is half the phase of a pure time delay 
of T seconds, giving rise to a rule of thumb tliat a DAC Is equivalent to a 
T/2 delay. This is always true for the phase characteristics. The gain 
characteristics, however, depend on input frequency, but this is an appro- 
priate approximation for •imall T. 
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Frequency analyzer- The frequency analyzer used was a Schlumbergcr 
model 1170. Its sine wave frequency generator provided the Input to the ADC. 
In addition, the analyzer used this input as well as the DAC's output to 
determine the gain and phase differences between the two. The appendix gives 
some Insight Into how the analyzer works in theory, but the Important thing 
from the system test point of view is the method by which the input sine wave 
is represented and Its interaction with the ADC sampling. Basically, the 
Schlumbergcr unit generates a sine wave in a discrete formulation. Two hun- 
dred discrete levels are used for every cycle. Thus the frequency generator 
effectively operates as though a continuous sine wave had been sampled, then 
as output through a zero-order hold device. 

The ADt; sampling the "discrete" sine wave is not synchronized in any way 
with the sampling rate of the frequency analyzer. This brings up two impor- 
tant considerations: the ratio of the sample-hold frequency of the analyzer 

to that of the ADC ..nd the observation that the ADC is, in effect, sampling a 
step input at each sample point. 

it was originally thought that the frequency ’■atio should be an Integer 
such as 10 or 20 merely for convenience and to ensure that the sampling fre- 
quency did not violate the Nyquist criterion of being less than twice that of 
the Ini'ut. However, the appendix shows that any spectrum generated by the 
ADC has no effect on the frequency analyzer operation when the frequency ratio 
divides evenly into 200. This leads to an operational requirement to change 
the input frequency when the sample rate of the computing system is varied. 

We let k • note that Tp^ * l/(200f), then f ■ k/(200T). 

This relationship specifies the input test frequency when k and T have been 
specified. The parameter k has the requirements that it divide Integrally 
Into 200 and that it be smaller than 100 to avoid violation of the Nyquist 
criterion. 

The fact that the sine w.ive is composed of 200 steps leads to a statisti- 
cal phase model of the interaction of the ADC and frequency generator. Fig- 
ure 5 Illustrates a typical A1>C sample of the input at some random point in 
the hold period, TpA. of the wave. Since the ADC sample may occur at any tlmi- 
in the period, a uniform distribution of pure time delays between 0 and Tp^ 
is expected. The mean delay is Tpy^/2. In terms of phase, the mean is 

*FA “ ”/200 radians 
- 0.9* 


kxperimf:ntal matrix 


Using the elements of the facl . y ar.d program as described here, the 
experimentation was designed to be ri-a in a manner akin to normal operational 
practices, that is, the values of T were representative, the software was 
that which is normally used, etc. The factors of the study were the step 
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FREQUENCY 

ANALYZER 

SIGNAL 


FA 




ADC SAMPLE 


COMPUTER CLOCK 
Figure 3.- Time delay due to frequency analyzer 
and ADC Interaction. 


size T, input frequency f, the particular digital filter (pure delay or 
first-order lag), and the presence or absence of the ADC pref liter. Table 1 
Illustrates these factors. 


TADLF. 1.- EXPF.KIME.NTAL FACTOPS 


Step size, T 

40, 5C msec 

Frequency, t 

1, 1.25, 2, 2.5 Hz 

Filter type 

e"®T. i/cs + 1) 

ADC prefilter 

in. Out 


The step size T was specified to be 40 and 50 msec in the faster calcu- 
l.itlon loop (80 and 100 ,isec in the slower loop). By letting It ■ 10 and 20, 
the frequency lest values were obtained. One might argue that k and not f 
is the factor: however, if the frequency had been selected first, the values 
of k would follow. It appears that either will suffice as the variable 
t actor . 

Care was taken with the first-order filter to ensure some ccnslstency 
within the test cases. The corner frequencies of the filters were set to 
1/10 the Nyqulst frequency or 

1/T - 2fl(l/2T)(0.1) 
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so chat 


T - IOT/tt 

While the value of T varies with T, it Is not a critical value in the .ense 
that results depend on it; therefore, it is not a factor in the experimental 
sense. The important aspect is that both classical first-order, low-pass 
filters and pure time delays were used to determine effects arising from time 
step size and time lags. 


RESULIS 


Data for the 32 cases resulting from the matrix were obtained and com- 
pared with theoretical or predicted results (table 2). In acquiring the phase 
data, it was immediately obvious that the results were statistical. The fre- 
quency analyzer technique yielded consistent run-to-run magnitude results, but 
the phase data varied considerably. Data for maximum and minimum phase values 
for each case were collected as a result; both are listed in tae table, except 
four cases for which data were not collected. The predicted values are the 
sums of the theoretical phases for the individual components in the hardware- 
software loop. The models for the DAC, digital filters, and prefilter were 
used with the values of T and f to obtain the prediction values. 

Comparing the predicted values with the experimental results shows that 
the experimental results follow fairly well with the deterministically pre- 
dicted values. From the models of the subsystems discussed here, only the 
frequency analyzer-ADC interaction and the possibility of time lags in the 
ADC-computer interaction are not included in the predicted values. The 
former is a phase-related item and the latter is a time lag item. This idea 
leads to a model for the difference or error between observed and calculated. 
Letting 


+ af 


where is the observed difference and Is the phasa offset that would 

represent phase-related effects. The product af is the portion of the error 
correlated with frequency or relative to a delay time since, for a pure time 
lag, the phase equals the product of frequency and the time lag. Figure 6 
shows phase errors plotted against frequency along with the linear regression 
fit to each of the two sets of data. Consider 

W = - 2-307f 


The constant indicates a phase offset of only 0.12°; this is not a 
and will be ignored as being within the error of data observation, 
the part correlated with frequency is Interpreted to result from a 
Tjo as 


large value 
However , 
time lag 
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TABLE 2.- THEORETICAL AND EXPERIMENTAL PHASE RESULTS 




















FREQUENCY, Hz 

Figure 6.- Linear regression fits to error data. 

-2.307f - 360“Tiof 
or 

T * 6*4 rnS6C 

Similarly, for • 0.09 - 1.6:3f, the lag is obtained as Tjq ~ ™sec 

Thus a lag of A to 6 msec apparently exists within the simulation loop, but 
no appreciable phase offset is apparent from the frequency analyzer-ADC 
interaction. 


CONCLUSIONS 


The goal of this study was to determine whether the VMS simulation com 
puter system could be verified using an experimental approach with a frequency 
analyzer. By carefully structuring the experiment and modeling the components 


the overall phase of the simulation loop Is both predicted and determined 
experimentally. The results Indicate that the technique can be used 
successfully. 

The differences between the predicted and measured values lend themselves 
to a statistical analysis which indicates that the computer hardware and soft- 
ware configuration imparts a time lt:g of approximately 5 msec to simulations 
on the facility. This value is relatively small since naive use of the 
facility can produce 8 msec by unintentional use of an ADC prefilter or as 
much as 40 to 50 msec with poor programming techniques. Such errors would be 
detected easily through the use of this technique. 
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APPENDIX 


EFFECTS OF ADC ON FREQUENCY ANALYZER 


The phase and gains produced by the frequency analyzer used in this test 
were derived on the basis that a sine wave input to a system will yield a sine 
wave output that may be represented as a sum of sine and cosine waves. Thus, 


Vq “ M sin(u)T + <>) 

» A cos (i)T + B sin u)T 


(Al) 


where 

M - (A^ + 


“ arctan A/B 


The frequency analyzer actually determines A and B first, then calcu- 
lates M and The procedure is directly from Fourier series analysis which 
yields 


and 


A 


B 


w 

n 


— 

n 



sin o)T 


cos U)T 


(A2) 


(A3) 


The frequency analyzer does not perform continuous integration but uses 
a simple quadrature formulation. The period of integration (summations) is 
one full cycle of the frequency. 

The ADC has the effect of creating copies of the signal input to it as 
described earlier in the text. For a sine wave input, its sampled version is 


V “ 


I 

T 


CD 



+ 




T 


The DAC effect on the spectra is to attempt to 
except n ” 0. However, the output contains traces 


filter out all components 
of the other components as 


sinTL + T + 
n-0 
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which may be written as 


Vo - I) [Ao + K + ^)] ^A4) 

n-0 '- 

This signal, Vg, Is the one applied to the frequency analyzer equations. 

Before performing the Fourier series integrations, It is useful to consider 
the particular frequencies used in the experiment. The computer frame time 
T was always chosen so that its frequency was a multiple of u). It was 
formulated as T - kn/lOOui. Thus the spectral locations are determined as 

. 2TTn /, . 200n\ /,c\ 

0) + — - 0)^1 + —^) (A5) 

The parameter k represents the number of computer sampling periods 
occurring for each period of the Input frequency u). Values of either 10 or 
20 used In this test make the frequencies of the spectra odd multiples of oj 
so that 

^(l + 2^) - (A6) 

where the 1^ are all odd Integers. When Vg from equation lAA) Is substi- 
tuted Into (A2), the fact that the Integral of odd functions over the Interval 
Is zero leaves only 

/ 2Tt/uj “ 

sin wT ^ aj^ sin igUii dt 
o n*0 

Using the fact that 1^ is odd leaves 

/ 2ti/u) 

3g sln^u)! - ag 
o 


Similarly, B 


Thus the frequency analyzer theoretically finds the amplitudes e; the 
sine and cosine of the primary frequency only. None of the repeated spectra 
are used to compromise the calculation of the frequency response of the s tern 
analyzed at frequency o. It must be stressed that this Is not a general 
result, but is true when the particular input sine wave frequencies and com- 
puter frame times are selected carefully. 
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